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1 Prologue 

In a talk where an elementary preamble to a theory of granular fluids was 
promoted pQ , one of us suggested a way to avoid observer dependence of some 
"thermal" like concepts: peculiar velocities should be read by reference to a 
frame translating with the speed of the centre of gravity (as is usual) , but also 
rotating again with an average speed dictated by the equation of balance of 
moment of momentum. Actually, because the Euler inertia tensor Y depends 
dramatically on the "affine" average speed, the "frame" is better chosen to 
deform affinely in accordance with equation of the balance of tensor moment of 
momentum. Thus one can avoid the introduction into the equations of the speed 
of rotation of the observer, as is done in expositions of extended thermodynamics 
and, more specifically for the present discussion, in [2] 

For a system of mass-points x^' (i — 1, 2...k) the relevant equations are (see 
(27) of d) 

fj& = / 

fi(K -BK-H) = M (1) 
Y = YB T + BY 
fi(H + BH + HB T ) = S 

where \i is total mass; x, place of centre of gravity; /, resultant of external 
forces /W; K is the tensor moment of momentum K = YB T ; B is the tensor 
rate of average affine displacement; M is the tensor moment of external forces 
Ylii x ^ ~ x ) ® f^i H 1S the Reynold's tensor evaluated on the peculiar speeds 
HH = M W « W ®« W with 

iW(r) = x(t) + G(r)«W(r), B = GG~ X \ 

S is the stirring tensor 

S = J2 G (s (i) ® G~ X f (i) + Gf _1 / (i) ® s {t) )G T . 

i 

* Department of Mathematics, University of Pisa 

t Department of Civil Engineering, University of Canterbury 



The extra complication attending with the splitting of the actual motion into 
the global translatory component x, the global affine component B{x^ — x), 
and the peculiar component Gs^ offers the advantage that quantities expressed 
in terms of the peculiar speeds «W are totally observer independent. Within the 
kinetic theory of gases, because disordered peculiar velocities are vastly larger 
than those ordinarily associated with the tensor B, the extra complication is 
not warranted; but that need not be so for granular gases and, more generally, 
for the class of 'kinetic continua'. 

Together with quantities mentioned above goes a kinetic energy tensor 

pW = ^fxx ® x + \^iBYB T + \-11E (2) 
and a corresponding kinetic energy theorem 

fjW= ~S + sym{x®f + BM). (3) 

If / and G~ Y M were constant and S were the total time-derivative of a "poten- 
tial" P, then a sort of principle of conservation of tensor energy would ensue: 

pW — -P + sym(x ® / + M) = constant. (4) 

A theory of "kinetic" continua may be based on the scheme imagined 
valid for material elements where p = p d(vol); p density. True, the densities 
of actions which correspond to /, M and S may not always be reducible to the 
Cauchy proposal (non-local action may be present). However, a first approach 
can be tried involving: external force per unit mass /; Cauchy's stress T; exter- 
nal tensor moment M and external stirring S both per unit mass; corresponding 
tensors of internal actions A and Z and third-order hyperstresses mm and as; so 
that relevant balance equations become 

p + p div x = 
Y -BY - YB T = 

px, = pf + div T (5) 
p(K - BK - H) = pM -A+ div m, 
p(H + BH + HB T ) = pS - Z + div s. 

Here S and Z must be symmetric and s must enjoy the property of minor left 
symmetry. 



Remark 1 . In classical fluid dynamics the gross interpretation of grad x or, 
better, of its symmetric and skew components is pervasive, random molecular 
motion notwithstanding. Still, to hide the effects of the latter motion outright 
within the thermodynamic maelstrom may, on occasion, curtail correct percep- 
tion of phenomena. A first coarse grasp of that recondite behaviour is offered 
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by the tensor H (a sort of tensor of deep ferment) through its "macroscopic" 
interpretation: write H in its canonical form highlighting eigenvalues (x^) 2 
and eigenvectors 

s=l 

(H being semidefinite positive its eigenvectors are non-negative, thus the choice 
of notation) and read it as follows: the population of grains (supposing that 
they all have the same mass) is spread among three tribes; within the s-th tribe 
the grains move along the line of h^ s \ their speeds (of magnitude x^) either 
pointing as the unit vector or as —h^ s ', the number of grains in each subtribe 
being equal. 

Alternatively one may imagine all grains to have (the same mass, say /i, and) 
the same speed intensity v but the fraction of those moving in the direction 
to be x /(Ex' s ' )• Under the latter circumstances v = p/fb is the number 
density of grains; if A is the mean free path, then (v/X)(p/p,) is the number 
density of collisions per unit time. 

Returning to the general case, we seem justified in calling v\~ 1 H 1 / 2 the 
collision density tensor. 

Remark 2 . The tensor H 1 / 2 has properties similar to those of the tensor 
introduced in an earlier paper of ours (for which, using different notation, the 
letter A was used, see sect.l of 0) in that H 1 / 2 n measures the "cross-over" 
rate through the plane of the normal n. That analogy apart, a fundamental 
difference remains because we have removed here the contribution to the peculiar 
speed of the affine motion and the Reynold's tensor of the earlier paper, call it 
H here, differs from H above if B does not vanish 

H = BYB T + H. 

2 Kinetic energy theorem. Balance of moment 
of momentum and sundry other balances 

The version for a continuous body B of the kinetic energy theorem Q, based 
here on density W per unit mass of the kinetic energy tensor 

W = -x <8> x + \bYB t + iff, 

is easily obtained by operating on Ipjlrrr with £®, on JSJ/v with B, summing 
term by term, taking the symmetric parts of each term, finally adding again term 
by term the last equation (J5J multiplied by \ and then, assuming smoothness, 
integrating over any subbody b of £>, by parts if need be, 
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f pW = f p[sym(i <g> / + BM) + is] + 

Jb Jb * 

f ■ 1 

+ / {sym[x ® Tn + B(mm)] + —xn} + 
Jdb 2 

[{-Z + sym[LT T + BA+bim*}}, (6) 

Jb 2 

where n is the unit normal vector to db, L = grad x,b = grad B, the exponent 
T indicates major transposition, and an exponent t to the third-order tensor m 
indicates minor right transposition. 

An ambiguity is left in the notation where third order tensors appear; rather 
than resolve the ambiguity by excess notation we leave it there, relying on the 
reader to sort it out easily. Just for once we declare that, in indicial notation 

(Bmn)ij = B ik m k jini, (bm*)^ = B ir<k m rjk . 

The sum of the first two integrals on the right hand side of © delivers the 
(tensor) power of external actions (respectively body forces, torques, stirring 
actions and boundary tractions, twists and ferment influx). The last integral 
must thus be interpreted as the tensor power of internal actions, of density 

- sym(^Z + LT T + BA+bmi 1 ); (7) 

hence the density of actual power is given by the scalar 

- [L -T + B ■ A T +b- {m t ) T + htr Z] (8) 

The fourth equation (JSJ , or rather its corollary obtained by taking the skew 
components of its two sides, though it exhausts the requirement of balance of 
vector moment of momentum, it does not here secure automatically the demand 
on stresses to make the internal power (JSJ observer- independent. Two observers 
on frames in relative motion read different values of L and B; the difference, in 
both, amounts to ew (e, Ricci's tensor; w, relative speed of rotation). Hence 
the condition 

skw T = skw A. (9) 

Oddly, it occurs sometimes that the constitutive choices for T and A are such 
that the stronger property 

T = -A T (10) 

applies. Then, as can easily be checked, even the the tensor power is 
observer-independent and reduces to 

-sym[^Z + (L — B)T T +tm*]. 
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Actually, in kindred investigations but where neither moments of momenta nor 
external torques are incorporated (K, M, ma all vanish in JSJ/y) one can dis- 
pense with a separate fashioning of the tensor A as that tensor would necessarily 
always coincide with pH again by our JSJl/y, now greatly reduced in content. 
Then, from QlOfl. the stronger identification obtains 

T = -pH. (11) 

Thus, through this constitutive law, a formal connection is enacted with pro- 
posals advanced in hypo elasticity, extended thermodynamics, etc., where the 
Cauchy stress is the main evolving function in an added balance equation. 

Another argument bears in favour of (|10fl . or, at least, reveals its deep gist. 
—A represents the density of internal equilibrated tensor torques; thus in the 
absence of twist influx due to subtler mechanisms it can be gauged in terms 
of T only as follows: imagine the material element as filling a minute sphere 
S t of radius e, obviously 'small' but not insignificant and thus imagine further 
— A(vol S e ) to be equal to the total over the surface of S e of the tensor moment 
of traction Tn (n, unit normal to S e ) 



-A{vol S e ) = en ® Tn. 
Js, 



If S e were the sphere of radius 1, so that f s n®n = ^ttI (I, the identity tensor), 
then 

^7re 3 A = -e 3 ( / n<8m)T T ; 



hence l(TUjl. 

Remarks . When (|1C)|> applies the tensor moment of inertia 
/ p{x ® x + K - BK - H) 

Jb 

is balanced by external actions only, i.e. by 



p(x (8 / + M) + / (x®Tn + mn), 

b Jdb 



for any subbody b. 



Another partial balance is often subsumed, at least as a constitutive prop- 
erty: the rate of change of total kinetic energy tensor 

( / pwy 

Jb 

is balanced by external tensor power only for all subbodies b. In other words the 
total tensor power of internal actions sums up to zero; its density J7J) vanishes. 
Then, necessarily 

Z = -2sym{LT T + BA+hm 1 }] (12) 
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Of course, such separate balances of powers need to be justified, if at all, 
by special physical circumstances. In particular relation H12|) applies at best in 
the absence of any dissipative effects or stirring effects from the macromotion 
and makes sense solely in conjunction with l|l(J|) . otherwise it would not be 
objective. Perchance both (|I0|l and l|12|l obtain only when written in terms of 
the conservative components of A and T alone. 

An argument similar to one called upon above to support l|10|) can be invoked 
in favour of the presence of some terms in i|I2|) . One needs only gauge also 
—Z(vol S e ) in terms of the virial of Tn over S e , the speed differential with 
respect to the centre being now (L — B)(en), 



hence the possible origin of the first two terms in (|I2|I . 

The expressions 0, have an important role in corollaries of the defini- 
tion of perfect internal constraints, i.e. of constraints such that the power of 
reactive internal actions vanishes for all virtual motions the constraints allow. 
For instance if the affine submotion at x is forced to coincide with that of the 
macromotion around x, i.e. if B = FF" 1 — L, then 



where the upper r is there to indicate the reactive contributions; (below an 
upper a indicates similarly active components). It follows that 



Under these circumstances the fourth equation JSJ) becomes irrelevant; at the 
same time the expression of T becomes 



If the definition of perfect constraint were to require the vanishing of the 



tensor power of reactions then Z would have to vanish, not only its trace. 

3 Boundary value problems; Constitutive laws 

The balance equations JSJ go along with appropriate conditions at the bound- 
aries which either render the constraints imposed there on x, B and H or em- 
body the local effects of the environment through the assignment of boundary 
traction Tn, twister mmn and stirrer sn. 

Actually, boundary conditions cannot be expected to mimic always the stan- 
dard model strictly. For instance, granularity and permeability of the restraining 





r rT r r 

T= -A , ™= 0, tr Z= 0. 



T = T+ [p(K -BK-H)- pM + A~div m] B =L- 



r 
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walls play sometimes a decisive role; their effects on the inner flow must be iden- 
tified and portrayed mathematically and that portrayal demands details on the 
'substructure' of the boundary. In any case the variety of continua for which 
the balance laws (0 are presumed to apply makes general statements unfeasible: 
loose granular matter is hardly entrained by a moving boundary or restrained 
by a stationary one, whereas no slip is allowed for viscous granular suspensions. 

Besides, the ingredient still missing is the set of constitutive laws for A and 
Z, T, m and s; each set characterises a member of the class of kinetic continua. 
Criteria of objectivity, thermodynamic compatibility, etc. restrict the choice of 
those laws, but we do not pursue the general issues here. Rather we pick a 
sufficiently comprehensive subclass, to encompass interesting even if disparate 
cases proposed in the literature and provide some explicit examples of flow. 

1. Standard stress is generated by deep ferment (like pressure in the kinetic 
theory of gases) 

T = -pH 

Such simple law applies for granular gases. However, most concepts and 
results of our analysis apply also, with some adjustments, to some vaguely 
similar settings, e.g. granular suspensions in a liquid. But then stress 
is influenced by viscous effects; additional terms enter the constitutive 
law for T with the involvement of L and B. Actually, requirements of 
objectivity rule out direct separate presence of these two tensors; they 
may enter only through the combinations 

D = sym L, sym B, L — B. 

The most elementary instance is when the dependence is additive and 
linear with some scalar coefficients of viscosity, say 771, 772, 773 respectively; 
when the medium is compressible further coefficients appear as factors of 
(tr L)I, (tr B)I, with I, identity tensor. Below the simplest occurrence is 
surmised: the usual linear dependence on D with viscosity rji as coefficient 
and an added linear dependence on L — B 

T = -pH + 2 m D + 2?/ 3 (L - B) 

thus viscosity may cause entrainment of the macromotion by deep twisting. 
A law of this type is suggested also in 0], see (5) of that reference. 

2. Again the simplest constitutive law is suggested for A which assures the 
validity of the mandatory condition JjJJl, but does not obey I|1U|I reflecting 
some of the preoccupations expressed just beyond equation (|12[) 

A = pH - 2 m (L - Bf . 

Insertion in equation ©jy shows that ferment does not influence K di- 
rectly, as the terms pH cancel out whereas twist connects macro and mi- 
cromotions. A reduced version, applicable to a reduced balance equation, 
is again suggested in 0], see (4) of that reference. 
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3. 



Twisting hyperstress is absent: 



ws. = 0. 

A linear dependence of am on grad B could be subsumed, by analogy with 

(6) of a. 

4. Relation 112fl applies provided dissipative contributions are first crossed 
out in T and A; in addition ferment is assumed to be, possibly, stifled 
by a sort of cross-over resistance (or collision loss) deemed to be propor- 
tional to pH and vice versa stimulated by the gross motion, the stimulus 
conjectured to be proportional to D 2 

Z — 2 sym[(L — B)pH\ + apH — 7Z? 2 ; a, 7 constant. 

5. The simplest rule applies for the stirring hyperstress: 

.$ = —f3 grad (pH); /?, a constant ferment transfer coefficient. 

Finally, the balance equations (J5J, in the local version, become 
dp 

— — h div(px) = 

OT 

QY 

— + (grad Y)x = BY + YB T 

OT 

I ' dx \ 

pi — +Lxj = pf + div [-pH + 2rnD + 2r) 3 (L-B)], 

( dB \ 
p[—+t>x + B 2 )Y = pM T + 2n 3 {L - B) T , (13) 



8t 

— + (grad H)x + LH + HL T 

OT 



P S + f3A {pH) - apH + ~fD 2 



The standard case is recovered, when 7 — 0, M — and S = 0, if one 
assumes that Y and H vanish, and cither B = L (and n\ then coincides with 
the usual viscosity) or B — skw L (and then n = n± + 773). If only Y were 
assumed to vanish, then —pH needs to be added to the viscous stress, H being 
a solution of the adjusted version of the last equation (|13l) . 

A preliminary reflection is appropriate: the fundamental law of moments 
fastens together main flow and twist forcing the constitutive laws to forge that 
link or else to pay the penalty of excluding skew components for T and A. On 
the contrary a possible direct connection of main flow with ferment, entrainment 
apart, is left to the hazards of the choice of constitutive laws, e.g. on the value 
of the constant 7 below. 

If alongside M and S also / vanishes, though 7 does not vanish, then one 
can seek a stationary solution of 113|) with constant density and L, B, Y, H 
constant tensors which must satisfy the following set of algebraic equations 

p(LH + HL T ) = -apH + 7D 2 , tr L = and 
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BY + YB T = 0, P B 2 Y = 2r] 3 (L - B) T . 

Besides, x must belong to the kernel of L: Lx = 0; so, as x — Lx, L 2 itself must 
vanish. Notice also that, as a consequence of the third condition, B 2 Y is equal 
to the symmetric tensor BYB T , hence 

pB 2 Y — 2773 sym(L — B), skw L — skw B. 

It is easy to verify that a solution exists where B, L and Y have only one 
non-null component (say B12, L12 and F 33 respectively) and H is determined 
consequently: with the choice above 



7 , L 2 
a— (! + - 2 



H13 = H23 — H33 — 0, Hu — 7^(1 H o")-^i2) 



•Hi 2 — — "J — o - ^12' ^22 — 



1 L 



4a 2 p 121 4ap 

and T is enhanced beyond the viscous contributions. 

Alternatively, and trivially, provided that a = 0, L and B may vanish alto- 
gether; x is then any constant vector and A any tensor field constant along the 
direction of x. 



4 Elementary flows 

All examples, except the last one, concern plane (£3 = 0, say) flows in an infinite 
channel: —00 < (1 < 00, < (2 < $, and with no external bulk influences: /, 
M, S vanish. 

Example 1. Granular gas with no loss in the bulk (a, /3/y, r\i vanish) nor 
losses on the walls. A simplest flow may be envisaged where stirring consists in 
a steady bounce wall-to-wall. If u is a constant vector along the first axis and v 
a similar vector but parallel to the second axis, one can choose initial conditions 
so that 

x = u, B = 0, H = v ® v. 

Granules jog up and down accross the channel with peculiar speed ±tt and, at 
the same time, move steadily down the channel. A bare image of the flow could 
be thus: at each point of the channel two clouds of granules meet, one with 
speed u + v and the other with speed u — v. 

The pressure exerted on the walls amounts to p\v\ 2 . 

Example 2. All conditions are as in the first example but in the presence 
of collision loss (a > 0). There exists a stationary solution where u is again 
a non-vanishing constant vector whereas v decays along the channel from the 
value «o at Ci = 0: 

x = u, B = 0, H = v (g)v e~ aCl 
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If u is the null vector (no flow down the channel), a solution independent of x\ 
exists where the bouncing between walls decays exponentially in time 



-ar 



H = vq ® vo e 

Example 3 aims to describe the effects of ferment loss due to collisions with 
the boundary by postulating that there the loss rate be proportional to ferment 
(7, a positive constant) 

— = - 7 if, at C2 = 0, C2 = S. 

OT 

A simple solution is found when there are no other losses as in Example 1, 
though ferment gradient affects the flow ((3 > 0). The solution involves the two 
constant vectors u, v as in the earlier examples, and two constants xo and 

x = u, 5 = 0, H = xo v®v e - Ca;i ~ 7T . 

All equations (|13(1 are trivially satisfied bar the last one which determines £ in 
terms of 7 

/?C 2 - MC + a - 7 = 0; 

an elementary discussion of subcases ensues, depending on whether the value of 
a falls within the interval 

or otherwise. 

Example 4. Plane Couette flow with constant imposed sliding velocity u 
in the upper plane £2 = S is one of the stationary flows hinted at (see end of 
Section 3): 

x = u((2/S); L = -u®c 2 , 


Ci, unit vectors along the axes. To the usual shear stress one must add 

7 M 2 / 2| U | 2 \ 7 \u\ 2 7 H 3 , 

4a (5 2 \ a 2 (5 2 / 4a <5 2 4a 2 <S 3 

If a and 7 go to zero, then only the constant component H±i remains arbitrary, 
all other components of H vanish. If only a vanishes, Hn is arbitrary again: 

rr lW\ 2 

and other components vanish. 

No connection can be expected between this example and the previous ones, 
as here no slip is allowed at the boundary. 
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5 Hints for progress 



So far we have striven to obtain evolution equations largely with the goal of 
portraying the behaviour either of granular gases (see, e.g., jo]) or of suspensions 
where the suspended granules are totally entrained by the surrounding viscous 
'solvent' although they provide the essential contribution to total inertia. Still, 
the balance laws © offer ground also for the study of the conduct of other 
continua; e.g., of hyperfluids, designed to model the evolution of more remote 
objects (see, e.g., 0). 

In this section we collect sundry remarks, handy when seeking links with such 
other pursuits. First, let us recall Remark 1 in Sect. 1; bearing the comments 
there in mind, we could say that our developments above (and any other based 
on J5J) alone) apply when the classification of granules within the three families 
quoted in that remark, though admittedly coarse, is nevertheless adequate. It 
seems unnecessary to go beyond it, if the bare aspects of global 'anisotropy' 
of the distribution of peculiar velocities are requisite: they can be evidenced 
already by the possibly different size of eigenvalues of H . 

Indeed, H offers a first appraisal of disorder in peculiar velocities; within 
an account in terms of H, maximum disorder is achieved when that tensor is 
spherical. Consider, for ease of display, circumstances where all granules have 
the same speed intensity v but random direction n; then 



where $(n) d{area) is the fraction of granules having speed in the immediate 
neighbourhood of vn and <S 2 is the unit sphere. When d is constant (thus equal 
to (47r) _1 ), H is spherical with the value 



hence we could introduce 




as an order tensor and repeat here developments formally identical with those 
available in the theory of nematic liquid crystals ,7,, in particular we could 
introduce a concept of tensorial temperance. 

The changes required so as to cover the case with generic speed intensity are 
now obvious: the distribution function "d{y) neds to be defined over J? 3 and to 
satisfy the normalisation conditions 




H = V/; 
3 ' 




then 
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and 



Q 



1 




(tr H) 



By analogy with the instance of perfect gases, some Authors have given special 
relevance to the canonical distribution 



thus evidencing a tensorial absolute temperature, as the inverse of the tensorial 
temperance O. 

Only further research will show if the field O is really decisevely better than 
the field of H in discussing real physical problems. Actually, there is a radical 
handicap in choosing the former: it is bound with the acceptance, without 
exception, of the canonical distribution, at least if the concept of temperature 
is itself not generalised, as vaguely muted in pQ, Sect. 5. 

Be that as it may, one exits here from the strictly mechanical province so as to 
admit thermodynamic (or, at least, thermodynamic-like) concepts, necessarily 
governed by the central axiom which expresses the balance of energy. Thus we 
close this section by proposing a formulation of that principle which may be 
appropriate for granular fluids. 

Actually we could perhaps dare to suggest that there should be, for strictly 
thermal phenomena in these continua, an inherent complexity parallel to the 
kinetic and dynamic one already imputed to them and consequently conjecture 
that it be possible to measure on each element at each instant along a process 
a density of thermal internal tensor energy, a third-order heat flux tensor, etc.; 
but we concede that such conjectures would be far-fetched here. We follow, 
rather, a middle course and, while suggesting a tensorial form of the principle of 
balance of energy to match the kinetic energy theorem (J3J, we take the deeper 
ferment to be isotropic and thus propose spherical tensors: to represent the 
thermal internal energy density and for the rate of heat generation; we also 
downgrade the third-order heat flux tensor to the form q ® I where q is the usual 
heat flux vector. Of course, as in the classical theory, e, A, q are linked with the 
'latent' molecular ferment, rather than the granular one argued about so far. 

In conclusion we postulate the validity, over any subbody, of a tensor balance 
equation modelled, formally, on the classical one 



/ sym(x ® Tn + B(imn) + —xn — q (£> n), (14) 
Jdb 2 

Along any process which is sufficiently regular to ensure the validity of the 
kinetic energy theorem (JBJ, equation (|14|> yields 
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and, because the choice of b among subbodies is arbitrary, the localisation ensues 
-pel = sym(LT T + BA+bm* - \-Z ■ -grad q + -pXI). (15) 

Taking the trace of both members a more common form of the energy principle 
is attained 

pe = L ■ T + B ■ A T +b- (m') T + \dr Z - div q + A. (16) 

As already remarked in Sect. 2 the first four terms in the right-hand sides of 
l|15|l and jTSjl assign the power (tensor or scalar) of internal actions; the last 
two measure heat loss or generation. 

Among the many consequences of (|16f) wc quote here, in conclusion, the 
following one: as mentioned before, when (|10l) applies, the first two addends in 
the right-hand side of l|16fl collapse into the product (L — B) ■ T. But 

L-B = G(G~ 1 F)F~ 1 , 

which suggest a conservative instance where e depends on G~ 1 F (as seems 
reasonable that it should) and 
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